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ABSTRACT 

We investigate the evolution of cosmological low mass (low virial temperature) 
objects and the formation of the first luminous objects. First, the 'cooling diagram' 
for low mass objects is shown. We assess the cooling rate taking into account the 
contribution of H2, which is not in chemical equilibrium generally, with a simple 
argument of time scales. The reaction rates and the cooling rate of H2 are taken from 
the recent results by Galli & Palla (1998). Using this cooling diagram, we also estimate 
the formation condition of luminous objects taking into account the supernova (SN) 
disruption of virialized clouds. We find that the mass of the first luminous object is 
several xlO^M©, because smaller objects may be disrupted by the SNe before they 
become luminous. Metal pollution of low mass (Ly-a) clouds also discussed. The 
resultant metallicity of the clouds is Z/Zq ~ 10"^. 

Subject headings: cosmology: theory — early universe — galaxies: formation 
— molecular processes — shock waves 
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1. Introduction 

Today, we have a great deal of observational data concerning the early universe. 
However, we have very little information about the era referred to as the 'dark ages'. 
Information regarding the era of recombination (with redshift z of about 10^) can be 
obtained by the observation of cosmic microwave background radiation. After the 
recombination era little information is accessible until z ~ 5, after that we can observe 
luminous objects such as galaxies and QSOs. On the other hand, the reionization of 
the intergalactic medium and the presence of heavy elements at high-z suggest that 
there are other population of luminous objects, which precedes normal galaxies. Thus, 
theoretical approach to reveal the formation mechanism of the such unseen luminous 
objects is very important. 

It is now widely accepted that luminous objects are formed from overdense regions 
in the early universe. These overdense regions collapse to form luminous objects, in 
case they fragment into many stellar size clouds and many massive stars are formed. In 
order to understand the way in which luminous objects are formed, physical processes 
of the clouds in various stages of evolution should be studied, individually. 

The formation process of luminous objects is roughly divided into three steps, 
formation of cold clouds by H and/or II2 line cooling, formation of the first generation 
stars in the cold clouds, and the star formation throughout the host clouds. These 
steps are disturbed by the feedback from the first stars. The first step has been 



investigated by many authors (e.g., Haiman, Thoul & Loeb 199(:, Ostriker & Gnedin 



T996| , [Tegmark et al. 1997| , [Gnedin fc Ostriker 199"7| , [Abel et al. T998D and it has been 



shown that the low mass clouds (virial temperature is several xlO^ K) become the 
earliest cooled dense clouds. The second step, however, is not investigated enough, 
although initial mass function and formation efficiency of the first generation stars 
in the clouds are very challenging and crucial problems. Many authors attacked this 
problem ( pVIatsuda, Sato fc Takeda 1969| , Hutchins 1976| , Carlberg 1981 , Palla, Salpeter) 



Stahler 1983| , Uehara et al. 199^ ) and they obtained various conclusions. However, 



now, the mass of the first generation stars are estimated through detailed investigation 
to be fairly large (Nakamura & Umemura 199£, Omukai fc Nishi 199^ ). For the third 



step, the feedback from the luminous objects on the other clouds has been studied 
by several authors ( [Haiman, Rees &: Loeb 1996 , Haiman, Rees & Loeb 1997, Ferrara 



1998| , paiman, Abel h Rees 1999| ). Haiman, Abel & Rees (1999) examined the 



build-up of the UV background in hierarchical models and its effects on star formation 
inside small halos that collapse prior to reionization. They stressed that early UV 
background below 13.6 eV suppresses the H2 abundance and there exists a negative 
feedback even before reionization. Moreover, the feedback from the formed stars on 
their own host cloud is more serious. The main feedback consists of two different 
processes, UV radiation from the stars and energy input by SNe. Through ionization 
of H (|Lin fc Murray 1992D and dissociation of H2 ( pilk 1977| , [Omukai h Nishi 19991) , 



- 3 - 



UV radiation have negative feedback on the further star formation in the host clouds. 
EspeciaUy, H2 is dissociated in such a large region that the whole of an ordinary low 
mass cloud is influenced by one 05 type star ( |Omukai &: Nishi 1999 ). The feedback 
from SNe on the host clouds is probably negative (e.g., [Mac Low &: Ferrara 199*9 ). 
SNe can disrupt the host clouds before they become luminous, because the explosion 
energy is comparable with the typical binding energy of host clouds. In this Letter, we 
investigate the evolution of low mass primordial clouds systematically, and assess the 
mass of the first luminous objects. 



2. Cooling diagram 

The formation of cold dense clouds, i.e., progenitors of luminous objects, is basically 
understood by the comparison between free-fall time and the cooling time. The 'cooling 
diagram' originally introduced by Rees & Ostriker (1977) and Silk (1977) shows the 
region where cooling time is shorter than the free-fall time on p - T plane, and vice 
versa. In this section, we present the cooling diagram on p - T plane including H2 
cooling. With this diagram, we can predict whether a cloud virialized at z = Zvir with 
virial temperature Tvir cools. 



2.1. H2 fraction with given virial temperature 

In order to estimate the cooling rate at T < lO'^K, we need the fraction of H2. The 
number fraction of H2 (hereafter denoted as 1/112) generally in equilibrium for 

T < 10^ K in the epoch of galaxy formation. The value of yu_2 at a time depends not 
only on p and T but also on the initial condition. Consequently, we cannot evaluate 
the cooling rate on the p - T plane without estimating non-equilibrium yn^. Tegmark 
et al. (1997) calculate yn^ numerically, however, their primordial yn^ is about two 
orders of magnitude over estimated because the destruction rate of by cosmic 
microwave background radiation at high- 2; is under estimated ( Galli &: Palla 1998| ). 
Since their primordial value ~ 10~^ is comparable to the necessary value to cool, 
their cooling criterion is not reliable generally. Here, using recent reaction rates and 
the cooling rate of H2 ( palli &: Palla 1998| ), we adopt a simplified and generalized 
method to estimate the cooling function of H2, differently from that of Tegmark et 
al. (1997). We introduce four important time scales, tdis, ifomn^cooi) and t^ec- They 
represent dissociation and formation time of H2, cooling time, and recombination time, 
respectively. Comparing these time scales, we assess the non-equilibrium fraction of H2 
with given virial temperature, and redshift. Below, our estimation is summarized (see 
Fig. |a). 

1. The case tdis < niin(tcooi, ^rec) (Region of "tdis fastest" in Fig. |l|a): H2 is in 
chemical equilibrium. In this case, = 1 where denotes the fraction of H2 in 
chemical equilibrium (solution of tform = ^di; 



12' 

lis)- 



- 4 - 



2. The case tree < min(tcooii ^dis) (Region of %ec fastest" in Fig. ||a): Ha is out 
of chemical equiUbrium, and H2 molecules are formed until the recombination process 
significantly reduces the electron fraction. As a result, yn^ is determined by the 
equation, tform = ^rec- Combined with the relation tform = ?/H2/j/Hi^dis( Susa et al. 
19981 ), is obtained as, = Uu^ 

{tree/ ^dis) • 



3. The case tcooi < iiiin(trec, ^dis) (Region of "tcooi fastest" in Fig. ||a): When the 
cooling time is the shortest of the three time scales, 2/H2 is determined by the equation 
iform = ^cooi- In Other words, increases until the system is cooled significantly. 



In case yH2 cooling dominates the other cooling processes, ?/h2 — y ^cooi/^dis- 
Otherwise, yu_2 is the solution of a quadratic equation. Here, t^^oi represents cooling 
time scale by H2 rovibrational transitions with y^. 

The electron fraction of a virialized cloud is assumed as ye = max(yg'^', Here 
yl*^^ is the fraction of cosmologically relic electrons calculated in Galli & Palla (1998). 
It equals to 3.02 x 10"^ for their standard model. The chemical equilibrium fraction of 
electrons is denoted yf^. With this electron fraction, we estimate the fraction of H2. 

In Fig. ||b, yH2 is plotted with given virial temperatures for four redshifts. For low 
redshift (z ^ 100) and high temperature (T ~ lO^K), H2 is in chemical equilibrium 
with given ionization degree. As the temperature drops, H2 gets out of equilibrium 
because the cooling time becomes shorter than the other time scales. Below 2000 K, 
recombination time scale is the shortest, and yu^ becomes relic value. 

For high redshift, the destruction of H~ {z ^ 100) and H^ {z ^ 200) by cosmic 
microwave background radiation reduces yn^ significantly. 



2.2. Comparison between free-fall time and cooling time 

We are able to assess the cooling rate with the given H2 fraction evaluated in the 
previous subsection. We compare the time scale of collapse (tg) with the cooling time 
(^cooi) which include the contribution from the H2 cooling. They are, 

iff = ( , tcool = l-5/»"^fcr.ir 

V32Gpvir/ ' ravirA(yH2 ) ^vir) ^vir) 

Here, p^ir = ISn'^^lpcr and rivir = 0,hPvir/mp, where pcr = 1.9x 10~^^/i^(l+Zvir)^ g cm~^. 
We adopt Q = 1, Qf, = 0.06 and h = 0.5 in this paper. 

Equating tg and tcool in eq. (|l]), we obtain the boundary between the cooled 
region during the collapse and the other region, which is drawn on the (1 + ^vir) - 
Tvir plane in Fig. |2|a. The objects virialized into the region denoted as tg < tcooi will 
be cooled by H2 during the gravitational collapse. In this case, collapsing cloud will 
be a mini-pancake, because the thermal pressure becomes negligible. We remark that 
the cooling region expands into Tvir ^ lO^K, which is different from classical cooling 
diagram such as the one in Rees &; Ostriker (1977). 
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We also compare the cooling time scale with the Hubble expansion time {H~^). 
The line icooi = is also drawn on Fig. |^a. The cooling region during the Hubble 
expansion time is slightly larger than the previous one, because Hubble expansion time 
is longer than the free-fall time. In this case, the collapse proceeds in semi-statically. 
As a result, the central region of the cloud will proceed to the run away collapse phase 
(iTsuribe k Inutsuka 1999| ). 



3. SNe and disruption of the bound objects 

As the collapse proceeds, small amount of the total gas is cooled to lOOK by H2. 



In those clouds, massive stars (10 — IOOMq) will be formed ( Nakamura fc Umemura 



1999 , Omukai &: Nishi 199§| ), eventually. After massive first generation stars form, 



evolution of the host clouds become slower because of strong regulation by UV 
radiation ( pmukai Sz Nishi 1998 ). Thus, next generation stars are hardly formed 



before the first generation stars die. Subsequent SNe might disrupt the gas binding 
before significant amount of total gas transferred into stars. Here, we derive the cloud 
disruption condition by SNe, with the assumption that the cloud is spherical and the 
density is constant, for the simplicity. 

We estimate the kinetic energy transferred from the SNe to gas. The velocity of 
expanding shock front from the center of a supernova remnant (SNR) is 

(2) 

where K = 1.53, E is the total thermal energy given by the SN, pi is the density of the 
cloud before the explosion, and t denotes the elapsed time since the explosion ( ppitzer 
1978). Integrating eq. (pi), we obtain the location of the shock front: 



'0.746(72 - 1) ke\ 



1/5 



R.{t) = t^l\ (3) 

The mass of the hot bubble is also obtained as msNR = |7ri?s (t)pi. The hot bubble 
keeps pushing the surrounding gas until the thermal energy is pumped off by the 
radiative cooling. Thus, the total momentum transferred from the SNR to the gas 
cloud is ptot = ''TiSNR(icooi)^^s(icooi)- Equating this momentum with the momentum of 
the whole cloud, we have the expanding velocity of the cloud: 

m'SNR(tcool) ^ 
Vtot = Vs (tcool) ■ (4) 

mtot 

If this velocity Vtot is smaller than the escape velocity (fesc) of the cloud, it will 
be still bounded. Otherwise, it disrupts. Now, we replace mtot in equation (^) with 
mj(Tvir, -Zvir)) which is the virialized mass of the cloud collapsed at Zvir with T^n- 
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The escape velocity from the cloud is directly related to the virial temperature. 
Consequently, we can draw the disruption boundary {vesc = "^tot) on the (1 + Zvir) - T^h 
plane. 

On the cooling diagram (Fig. |2|a), the boundaries {vesc = Vtot) are superimposed 
for two cases. These lines are obtained with the assumption that the input thermal 
energy from the SNe is 10^^ erg and 10^^ erg, respectively. The input energy almost 
reflects the number of the SNe. The values 10^^ erg and 10^^ erg represent the case of 
single SN and 10 SNe, respectively. The former should corresponds to the clouds in 
tg < tcooi < H~^, because they will have a runaway collapsing central core. The core 
evolves much faster than the envelope and will be a massive star, probably, followed 
by a single SN. The latter case represents the clouds in tg > tcooi- They will have 
a shocked pancake, and the cooled region will fragment into stars. That's why they 
should have multiple SNe. 

However, we should note that the disruption criteria by SNe strongly depend on the 
geometry of objects ( [Mac Low &: Ferrara 1999| , piardi et al. 1999 and references there 
in). In the case that cooling is efficient {tg > tcooi)) a cloud evolves dynamically and 
becomes complicated shape, which is probably flattened. The geometric effect makes 
the momentum transfer from the SNR to the surrounding gas less efficient than our 
evaluation. On the other hand, if cooling is not efficient {tg < tcooi)) a cloud becomes 
fairly spherical and has a centrally condensed density profile. In this case, the effect 
of SNe may become stronger than the above estimate (e.g., Morgan & Lake 1989). 
Moreover, if the duration of the multiple SNe is longer than or comparable with the 
evolution time scale of a SNR, disruption criteria is not evaluated only with the total 
energy of multiple SNe (e.g., Ciardi & Ferrara 1997). Thus, our estimate shows the 
qualitative tendency, so that detailed calculation for a individual cloud is necessary to 
derive the SNe effects accurately. 



4. Evolution of low mass objects and mass of the first luminous objects 

According to the argument in the previous section, the clouds in the region 
tg > tcooi will experience multiple SNe. Assuming that the total energy of multiple SNe 
as 10^^ erg, the survived region is the dark shaded upper region of Fig. |2|b (denoted 
as "LO").0 As shown in Fig. ^3 survived clouds are fairly massive {T^n ^ 10^ K) 
and they evolve into luminous objects through following processes (iNishi et al.l^98|) : 
(1) By pancake collapse of an overdense region or collision between subclouds in a 
potential well, a quasi-plane shock forms (e.g., Susa, Uehara fc Nishi 1996| ). (2) If 



the shock-heated temperature is higher than ~ 10 K, the post-shock gas is ionized 



3 



Of course there exists some ambiguity in the total energy, but Fig. Qb dose not change by this ambiguity 



quaUtatively. 
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and cooled efficiently by H line cooling. After it is cooled bellow 10^ K, H2 is formed 



fairly efficiently and it is cooled to several hundred K by H2 line cooling (e.g., Shapiro 



Kang 1987, |Susa et al. 1998 ). (3) The shock-compressed layer fragments into 



cylindrical clouds when tdyn ~ ifrag ( Vamada &: Nishi 1998 , Uehara fc Nishi 1999 ). 



(4) The cylindrical cloud collapses dynamically and fragments into cloud cores when 



^dyn ~ ^frag ( Uehara et al. 1996 , Nakamura &: Umemura 1999 ). (5) Primordial stars 



form in cloud cores ( Omukai &: Nishi 1998 ). (6) Since the gravitational potential of 



the cloud is deep enough, subsequent SNe cannot disrupt the cloud. Star formation 
regulation by UV radiation is also weak because of highly flattened configuration of 
the host cloud. (7) Next generation stars can form efficiently and the cloud evolves 
into a luminous object. 

On the other hand, the clouds in the region < tcooi < will experience a 
single SN. As a result, the survived region is bounded by the line denoted as 10^^ 
erg (denoted as "Ly-a" of Fig. ^) and they evolve into luminous objects if they are 
isolated. However, the evolution time scale of these clouds is rather long because of 
large disturbance by the SN and they are reionized at low-z. Q After reionization, 
they may be observed as Ly-a clouds. Since the baryonic mass of these clouds are 
several xW^Mq and the ejected metal mass by a SN is several Mq, their metallicity 
is estimated to be ~ 10^'^ Zq. The observation of QSO absorption line systems imply 



the similar metallicity for Ly-a clouds (Cowie et al. 1995; Songaila fc Cowie 1996 



bong aila 19971 ; |C owie fc Songaila 1998| ). They are typically the la objects and collapse 
at z ~ 10. Our estimate basically agrees with the more detailed calculation in Ciardi 
k Ferrara (1997). 

In the unshaded lower right region (Tcmb > T^ir) and the lightly shaded region 
"NC" of Fig. clouds are diffuse and do not become luminous, because radiative 
cooling is not efficient. In the shaded region "IG", SNe destroy the binding of host 
objects, followed by the diffusion of heavy elements into the surrounding medium. 

Therefore, the ffist luminous objects are probably formed in the region "LO" and 
their mass is estimated to be several xIO'^Mq, if we consider the 2 ~ 3cj objects. 
Formation epoch of the first luminous objects is 2; ~ 30 (considering the 3cr objects) 
or z 20 (considering the 2a objects). This estimated mass is larger than the one 
obtained by Tegmark et al. (1997), because small clouds may be blown up by their 
own SNe. 
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Fig. 1. — (a) The Ue - T plane is divided into three region. In the region with higher temperature, 
idis is the fastest and yii^ = y'^^. In the region with lower temperature, tree is the fastest and 

= yn^l^- Between these two regions, where tcool is the fastest, yji^ 2± VuiV'td^' (^) Virial 
temperature v.s. H2 fraction is plotted for given epoch of collapse (z). 
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Fig. 2. — (a) Cooling diagram is plotted on the (1 + z^ii) - TVir plane. The thick solid lines divide 
the plane into three regions. In the region denoted as tfi > icoob cooling proceeds faster than 
the collapse. The narrow region denoted as ts < tcooi < H~^, the cloud cools faster than the 
Hubble expansion, but cannot catch up with the free-fall. The clouds virialized into the lower 
region {H^^ < tcool), cannot cool and will not be the luminous objects. The lower right region 
(TciMB > T^it) is the forbidden region of cooling, due to the Compton heating. Two dashed lines 
denote the boundary below which virialized clouds are destroyed by SN explosions. The lines 
correspond to the cases that the total energy of SNe is 10^^ erg and 10^^ erg. (b) (1 + Zyir) - ?Vir 
plane is divided into five regions. In the unshaded lower right region (Tcmb > T^ir) and the lightly 
shaded region denoted as "NC" , clouds are diffuse and do not become luminous, because radiative 
cooling is not efficient. In the shaded region denoted as "IG", clouds will be disrupted by the 
SNe. Clouds in the region "Ly-a" will not be destroyed by the SN. The dark shaded upper region 
(denoted as "LO"), represents the clouds which could be luminous objects. The thick solid three 
lines marked as la, 2a, and Sa represent the location that the density perturbations in standard 
CDM cosmology should virialize. The mass of the clouds are constant along the thin solid four 



